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An  upper  bound  on  the  height  of  the  peak 
sidelobc  of  the  random  array  is  derived  using 
the  theory  of  level  crossings  of  a random  pro- 
cess.  This  bound  is  obtained  by  calculating 
the  expected  nun^er  of  uperossings  of  a given 
level  by  the  array  power  pattern.  The  bound  is 
compared  with  computer  simulation  results. 

I . INTRODUCTION 

The  linear  random  array  has  been  studied  by 
several  authors,  e.g.  [l]-[4].  In  this  type  of 
array  all  of  the  antenna  elements  lie  on  a line, 
and  their  positions  are  chosen  randomly  and  in- 
dependently from  some  probability  distribution. 

So  any  array  with  specified  element  locations 
designed  in  this  way  is  a particular  realization 
from  the  set  of  possible  arrays  determined  by 
the  given  probability  distribution.  By  choosing 
positions  in  this  way,  fewer  elements  may  be 
used  than  would  be  necessary  in  an  array  with 
uniformly  spaced  elements,  and  grating  lobes, 
which  are  sidelobes  of  height  equal  to  the  height 
of  the  main  lobe,  will  be  eliminated  111.  Be- 
cause of  this  property  a random  array  could  be 
useful  in  an  application  in  which  a few  elements 
must  be  spread  over  a large  aperture.  Since  the 
element  locations  are  chosen  randomly,  the  power 
pattern  of  the  array  is  a sample  function  of  a 
random  process  which  is  a function  of  angle.  It 
Is  important  to  be  able  to  characterize  the 
Plight  of  the  peak  (i.e.,  highest)  sidelobe  of 
jlhe  power  pattern.  A complete  characterization 
of  the  height  of  the  peak  sidelobe  would  be 
given  by  its  probability  density  function  (pdf) . 

Approximations  to  this  probability  density 
function  Juve  been  derived  (l]-{4]  by  consider- 
ing samples  of  the  power  pattern  at  a number  of 
points.  Different  results  will  be  obtained  in 
this  paper,  where  an  upper  bound  on  the  height 
of  the  peak  sidelobe  will  be  derived  by  using 
the  theory  of  level  crossinas  of  a random  pro- 
cess (discussed,  for  example,  in  [5]-[7]). 

Since  it  is  important  in  antenna  array  design  to 
make  the  sidelobes  as  low  as  possible,  a proba- 
bilistic upper  bound  on  the  peak  sidelobe  height 
can  be  of  considerable  significance. 


is  the  same,  then  S^O) 


exp{ jky^sin0)  is  the 


8(0^, 6)  = 


I exp[ jky. (sin0  - sin0-)]. 
i=l  ^ 


The  power  pattern  18(0^, 8) | 


of  the  array  is  the 


power  received  from  a sinusoidal  source  at  angle 
0 when  the  array  is  weighted  to  look  in  the  dir- 
ection 0Q. 

Since  the  element  locations  are  random 


quantities,  B{0^,6)  and  |b(0q,0)| 


are  random 


processes. 


The  special  case  in  which  0q*O  will 


be  considered,  although  an  extension  is  possi- 
ble to  the  case  in  which  0q?^O,  Let  u=sin0, 

X(u)  = B(O,0)  , and  let  Xj^(u)  and  ^^(u)  be  the 

real  and  imaginary  components  of  X(u).  Note 
that  u lies  in  the  interval  [-1,1].  From  (1), 

X(u)  = X^(u)  + jX2(u) 


II.  SIDELOBE  CHARACTERISTICS  OF  THE  RANDOM  ARRAY 

Consider  a linear  antenna  array  with  N ele- 
ments distributed  over  an  aperture  of  length  L. 
Let  y^»y2* • • • ^ randomly  chosen  loca- 
tions of  the  array  elements  with  respect  to  an 
origin  arbitrarily  picked  at  the  midpoint  of  the 
aperture.  Suppose  the  wavelength  of  incident 
radiation  is  X,  for  which  the  wave  number  is 
k=2iT/A.  If  the  power  received  at  each  element 


complex  amplitude  of  the  signal  received  at  the 

i^*^  element  from  a sinusoidal  point  source  at 
angle  0 from  array  broadside.  The  complex  weight 

needed  to  bring  the  i^^  element  signal  into  phase 
with  the  signal  at  the  origin,  so  that  the  sig- 
nals may  be  coherently  summed,  is  W^(0)  = 

exp(-jky^sin0) . The  complex  beam  pattern  B(0q,0) 

is  defined  to  be  the  complex  amplitude  of  the 
signal  received  by  the  array  from  a sinusoidal 
source  in  the  direction  0 with  element  weights 
fixed  at  W.  (0,.).  Thus 


This  upper  bound  will  be  compared  with  com- 
puter siraulations.  It  will  be  shown  that  the 
derived  result  yields  a reasonably  tight  upper 
bound  on  the  peak  sidelobe  height. 


I cos(ky  u)  + j 
i-1  ^ 


t sin(ky .u) , (2) 

i«l  ^ 


bound  on  the  peak  sidelobe  height.  and  the  power  pattern  is  |x(u)  | • Since  the 

beam  x>attern  X(u)  has  the  syimetry  property 
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X(-u)  = X*(u),  X(u)  in  the  interval  (0,1)  deter- 
mines X(u)  in  the  interval  (-1,0).  So  it  is  suf- 
ficient to  consider  the  beam  pattern  X(u)  only 
over  the  interval  (0,1). 

The  random  process  X(u)  can  be  partially 
described  in  terms  of  the  means,  the  autocovari- 
ance functions,  and  the  crosscovariance  function 
of  Xj(u)  and  x^iu).  The  special  case  in  which 

the  element  positions  are  independent  and  uni- 
formly distributed  over  the  interval  (-L/2,L/2) 
will  t>e  considered.  The  mean  value  of  the  pro- 
cess X^(u)  is  (4) 

E{Xj^(u)}  = N sinc(uL/X)  (3) 

for  ue(0,l),  where  sinc(uL/A)  = . The 

TTUL/A 

mean  value  of  is  e{x2(u)}  = 0- 

It  can  be  shown  that  the  crosscovariance 
function  K (u, ,u-)  is  zero  for  all 

U2,U2c[0,l].  In  the  sidelobe  region  of  the  beam 

pattern,  where  u is  much  larger  than  X/L,  the 
autocovariance  functions  of  X^Cu)  and  X2(u)  are 

approximately  equal  and  are  given  by 

\ (v)  = (v)  = ^ sinc(vL/X),  (4) 

1 ^2 

where  v=U2-U2,  The  mean  of  the  process  (u)  can 

be  neglected  in  the  sidelobe  region  if  it  is 
small  in  comparison  with  the  standard  deviation 

/n/2  of  the  processes.  The  mean  of  X^Cu)  will 

first  be  neglected,  and  later  its  effect  on  the 
peak  sidelobe  will  be  considered.  Therefore, 
here  it  will  be  assumed  that  the  processes  X2(u) 

and  X2(u)  are  wide-sense  stationary  in  the  side- 
lobe region  of  the  beam  pattern.  The  processes 
X2(u)  and  X2(u),  for  a given  u,  are  each  equal 

to  a sum  of  N independent,  identically  distrib- 
uted random  variables.  A multivariate  form  of 
the  central  limit  theorem  [8]  can  be  used  to 
show  that,  in  the  sidelobe  region,  X2(u)  and 

X2(u)  can  be  approximated  to  be  normal  random 

processes  with  zero  means  and  autocovariance 
functions  given  by  (4) ; henceforth  they  will  be 
treated  as  if  they  were  normal  processes.  Since 
X^Cu)  and  X2(u)  are  assumed  to  be  normal  and  are 

uncorrelated,  they  must  be  independent,  and  since 
each  is  wide-sense  stationary,  they  must  both  be 
stationary. 

The  peak  sidelobe  height  is  the  highest 
local  maximum,  other  than  the  maximum  at  u-0 
which  is  the  main  lobe,  of  the  power  pattern 

|x(u)|  « xjcu)  + X2(u).  The  expected  number 

of  uperossings  of  a level  z by  the  square  root 


t 


V2  2 

X^Cu)  + X2(u)  of  the  power  pattern,  in 

any  interval,  can  now  be  calculated.  This  will 
yield  a bound  on  the  probability  of  K(u)  exceed- 
ing z in  the  sidelobe  region.  It  can  be  shown 
that  the  expected  number  E{u^)  of  uperossings  in 

an  interval  of  length  % of  the  level  z by  any 
stationary  random  process  A(u)  is  given  by  (7) 


E{U  } f f |a  |f..'(Zfa  )da  (5) 

2 jo  ^ 

where  f^i  (a, a ) is  the  joint  pdf  of  the  pro- 

I 

cess  A(u)  and  its  derivative  A (u) . (It  will 
be  assumed  here  that  A(u}  has  continuous  sample 
functions.  The  power  pattern  of  the  random  ar- 
ray has  this  property) . This  expectation  value 
has  been  calculated  for  the  case  in  which  A(u) 
is  normal  [5]  and  the  case  in  which  A(u)  is  the 
envelope  of  a normal  process  (6).  Methods  used 
by  Rice  [5]  will  be  used  here  to  determine 

for  R(u).  The  joint  pdf  of  X.  (u)  and  X.(u),  and 

• t ^ 

their  derivatives'  X^Cu)  and  X2(u),  at  any  point 
u,  is 


I I 

' ’ (x, ,X-,X, ,x,)  = 


Xj,X2,Xj^,X2  1'  2'  1'  2 


exp 


(2n) 


r , 2 2 
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■ 2^ ' 2 n 

2>  X ^ 
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( 0 
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. 

(6) 


where  X =N/2  is  the  variance  of  X (u)  and  (u) , 
2 * 

and  * (N/6)  (ttL/X)  is  the  variance  of  X^  (u) 

and  X^tu).  These  variances  can  be  calculated 

using  the  technique  described  in  [6],  for  ex- 
ample. By  using  the  change  of  variables  = 

r cos6,  X2  = r sin6,  the  joint  pdf  of  R(u)  = 

^X^ (u)  + X2 (u) , 0(u)  = tan  ^ [X^ (u) /X2 (u) ] , and 

their  derivatives  R (u)  and  6 (u) , cm  be  found. 
Integrating  over  all  possible  values  of  6 and  9* 
yields  the  joint  pdf  of  R(u)  and  R* (u) s 
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•^^X^  Xq 


(7) 


Using  equation  (5)  for  the  expected  number  of 
uperossings  of  z by  R(u),  and  setting  il  = 
(1-X/L) , the  length  of  the  sidelobe  region, 
gives 

E{U^)  = (1-X/L) 


■ Stetta 
MeUpn 


= (l-X/U  ^ t i ^ (8) 
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Equation  (8)  is  the  main  result  of  this  section, 
and  will  be  used  to  set  an  upper  bound  on  the 
height  of  the  peak  sidelobe. 


e(0  } = [ du[  |r  (u) |f  > (z,r' (u) ;u)dr  , (9) 

* Jx/b  >0 


To  use  (8)  to  study  the  peak  sidelobe  of 
R(u),  the  relationship  between  level  crossings 
of  a random  process  and  local  maxima  of  that 
process  must  be  considered.  There  are  at  least 
as  many  local  maxima  above  z as  there  are  up- 
crossings  of  z;  thus  the  probcibility  of  an  up- 
crossing  of  z is  equal  to  the  probability  of 
having  at  least  one  local  maximum  above  z. 
Following  methods  used  by  Rice  [5],  it  can  be 
shown  that  the  expected  number  of  local  maxima 
above  z approaches  the  expected  number  of  up- 
crossings  of  z for  high  levels  of  z.  This 
justifies  the  use  of  level  crossings  to  study 
the  local  maxima  of  R(u}. 

For  the  case  of  R(u)  the  expected  value 
E{U^}  of  the  number  of  uperossings  in  the  side- 
lobe region  is  an  upper  bound  on  the  probabil- 
ity of  at  least  one  uperossing  of  z in 

that  region.  The  number  of  uperossings  of  a 
level  z can  assume  only  non-negative  integer 
values.  Letting  p^(z)  be  the  probability  of 

exactly  i uperossings  of  z,  it  is  clear  that 
e{u  } = r ip.(z),  and  P{u  } = Z p. (z) . 

Z . ..  L Z 1 

1=0  1=1 

Therefore  ^ e{u^},  and  upper 

bound  on  It  is  reasonable  to  expect 

that  this  bound  becomes  tight  for  small  values 
of  E{u^},  since  in  this  case  the  probability 

of  obtaining  more  than  one  uperossing  in  the 
sidelobe  region  is  negligible  compared  to  the 
probability  of  obtaining  zero  or  one  uperossing. 

The  case  in  which  the  mean  e{X^(u)}  of 
(u)  can  be  neglected  has  been  considered. 

Since  e{x^(u)}  = N sinc(uL/l)  and  the  standard 
deviation  of  Xj^(u)  is  approximately  /n/2,  we 

see  that  as  N increases  while  L/X  remains  con- 
stant, the  mean  will  increase  relative  to  the 
standard  deviation,  and  wil]  make  a large  con- 
tribution to  the  beam  pattern.  The  case  in 
which  e{x^(u)}  cannot  bo  neglected  will  now  be 

treated.  The  process  X^(u)  is  no  longer  sta- 
tionary, and  the  expected  number  of  uperossings 
of  z by  R(u)  is  not  given  by  (5) . 

The  expected  number  of  uperossings  by  an 
arbitrary,  nonstationary  random  process  can  be 
calculated  by  integrating  over  the  interval  in 
question.  Specifically  the  expected  number 
E{U^}  of  uperossings  in  the  interval  [A/L,l]  of 

the  level  z by  the  square  root  R(u)  of  the  power 
pattern  is  given  by  ([6],  Chapter  13) 


where  fj^*(*’tu),r  (u);u)  is  the  joint  probabil- 
ity density  function  of  R(u)  and  K (u) , written 
as  a function  of  ue[X/L,l).  (Note  that  u is  not 

a random  variable) . The  expected  value  E{U  } is 

z 

not  obtainable  in  closed  form,  and  must  be  cal- 
culated by  numerical  integration.  Computed  re- 
sults based  on  (9)  will  be  presented  in  graphi- 
cal form  in  the  next  section. 


III.  COMPARISON  WITH  COMPUTER  SIMUIATION  RESULTS 

The  bounds  on  the  peak  sidelobe  height 
given  by  equations  (8)  and  (9)  are  presented 
graphically  in  Figures  1 and  2,  where  the  solid 
lines  represent  the  bound  given  by  (9)  and  the 
dashed  lines  those  given  by  (8) . The  value  of 
2 2 

the  normalized  level  z /N  versus  the  normalized 
length  L/X  of  the  array  is  plotted  for  different 
numbers  of  elements  N,  and  the  indicated  values 

2 2 

of  E{U  }.  The  normalized  level  z /N  is  used 
z 

because  from  it  is  obtained  the  height  of  the 
peak  sidelobe  of  the  power  pattern  relative  to 
the  height  of  the  main  beam.  In  Figure  1,  the 
value  of  E{U^}  is  0.1,  and  in  Figure  2,  e{u^} 

is  0.3.  From  Section  II,  we  know  that  the 
probability  of  at  least  one  uperossing  of 

z is  less  than  or  equal  to  E{u  },  Thus  for  a 

z 

given  N and  L/X,  the  probability  P{u  } of  the 

z 

normalized  height  of  the  peak  sidelobe  of  the 

2 2 

power  pattern  exceeding  the  value  z /N  given 
by  the  curve  of  Figure  1 or  2 is  less  than  or 
equal  to  E{U^}.  This  is  equivalent  to  saying 

that  the  normalized  height  exceeded  by  the  peak 
sidelobe  with  probability  E(U  } is  less  than  or 

* 2 2 

equal  to  the  normalized  height  z /N  given  by 
the  curve  of  Figure  1 or  2.  It  can  be  seen  that 
there  are  differences  between  the  calculations 
based  on  the  assumption  that  the  beam  pattern  is 
stationary  with  zero  mean  and  the  more  precise 
calculation  which  considers  the  nonconstant  mean 
of  the  beam  pattern.  The  peak  sidelobe  upper 
bound  increases  due  to  the  mean  term,  particu- 
larly where  the  density  N/(L/X)  of  array  ele- 
ments per  array  length  is  high.  It  is  reason- 
able to  expect  this  difference  to  increase  when 
N increases  or  L/X  decreases,  since  in  these 
cases  the  mean  term  becomes  larger  with  respect 
to  the  standard  deviation  of  the  beam  pattern 
over  the  sidelobe  region. 

Figures  1 and  2 also  contain  the  results  of 
computer  simulations  that  were  done  to  study  the 
validity  of  the  upper  bounds.  For  given  values 
of  N and  L/X,  a set  of  random  arrays  was  formed 
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by  gcnGratin9  independent  random  numbers  to  rep- 
resent element  positions  uniformly  distributed 
over  (-L/2,L/2],  The  power  pattern  of  each  ar- 
ray was  calculated,  and  the  peak  sidclobc  was 
found.  For  each  value  of  N and  L/X,  an  empiri- 
cal cumulative  distribution  function  of  the 
peak  sidelobe  value  was  formed  from  this  set  of 
arrays.  In  Figure  1,  the  points  represent  Uie 
normalized  level  which  exceeded  the  peak  side- 
lobe  height  of  90%  of  the  arrays,  and  in  Figure 
2,  the  points  represent  the  level  which  exceed- 
ed the  peak  sidelobe  height  of  70%  of  the  ar- 
rays. For  each  value  of  N and  L/X  shown,  100 
arrays  were  generated,  except  for  N=20,  L/X=200, 
for  which  50  arrays  were  generated.  However, 
for  each  N and  L/X,  the  points  in  Figures  1 and 
2 are  taken  from  the  same  empirical  distribution 
function,  and  do  not  represent  the  results  of 
inde}>cndent  experiments.  It  can  be  seen  that  in 
all  cases  the  simulations  fell  below  the  curve 
given  by  the  more  precise  upper  bound  obtained 
from  (9) , and  for  lower  densities  of  elements 
N/(L/X),  the  experimental  points  fell  below  the 
upper  bound  calculated  from  (8) . 

II  is  likely  that  the  density  of  elements 
will  not  be  high  for  practical  random  array 
designs.  For  an  antenna  array  with  a density 
of  N/(L/X)=2  or  greater,  a periodic  array  can 
be  designed  with  no  grating  lobes  and  low  side- 
lobes,  and  for  arrays  with  such  densities, 
mutual  coupling  between  elements  changes  the 
beam  pattern,  making  equation  (1)  invalid.  Us- 
ing random  element  positions  would  be  most  use- 
ful in  applications  where  a low  density  of  ele- 
ments must  be  used,  since  in  that  case  the  ef- 
fect of  mutual  coupling  is  small,  and  the  use  of 
equally  spaced  elements  could  lead  to  a large 
number  of  grating  lobes.  In  these  cases,  the 
upper  bound  given  by  (8)  could  be  used  to  deter- 
mine the  peak  sidelobe  height.  This  is  fortunate 
since  (8)  is  simple  to  compute.  Using  (9)  to 
determine  the  upper  bound  would  only  be  necessary 
in  the  loss  practical  case  of  a high  element  den- 
sity. 

In  Section  II,  the  approximate  covariance 
function  (1)  was  derived  based  on  the  assumption 
that  the  region  near  the  main  beam  was  excluded 
from  the  silelobe  region.  The  region  (X/L,l] 
has  been  considered  to  be  the  sidelobe  region  for 
the  calculations  and  simulations  presented  in 
this  section,  but  the  results  might  have  been 
more  valid  if  a smaller  sidelobe  region,  such  as 
(3X/L,1],  had  been  used.  The  region  [X/L,l]  was 
studied  because  for  practical  array  design,  it  is 
necessary  to  determine  the  peak  sido?obe  over  the 
entire  region,  not  only  over  the  region  in  which 
.the  Uieory  yields  a very  accurate  result.  Since 


the  region  near  the  main  beam  exhibits  some  of 
the  properties  of  deterministic  arrays  (4),  an 
alternative  method  of  studying  the  peak  sidelobe, 
which  will  not  be  considered  here,  might  be  to 
use  the  theory  to  study  the  restricted  region  in 
which  the  approximations  arc  very  accurate,  and 
treat  the  region  near  the  main  beam  as  if  it 
were  deterministic. 

Figures  3 and  4 compare  the  results  derived 
in  this  paper  with  previous  results  concerning 
the  peak  sidciobe  of  the  random  array.  Steinberg 
[3],  [4]  used  a technique  in  which  the  beam  pat- 
tern was  sampled  at  a number  of  points,  an  esti- 
mate of  the  peak  sidelobe  height  was  formed,  and 
a correction  term  based  on  interpolation  was 
added  to  the  estimate.  The  solid  and  dashed 
lines  and  simulation  points  in  Figures  3 and  4 
were  taken  from  Figures  1 and  2,  respectively, 
and  the  dotted  line  was  calculated  using 
Steinberg’s  technique.  It  can  be  seen  that 
Steinberg’s  estimate  is  equal  to  or  greater  than 
the  upper  bound  derived  in  this  paper,  and  lies 
above  the  simulations  reported  in  this  paper, 
for  the  parameter  values  that  were  considered. 
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